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Abstract: We study the viability of a nonlocal dispersal strategy in 
a reaction-diffusion system with a fractional Laplacian operator. We 
show that there are circumstances - namely, a precise condition on the 
distribution of the resource - under which the introduction of a new 
nonlocal dispersal behavior is favored with respect to the local dispersal 
behavior of the resident population. 

In particular, we consider the linearization of a biological system that 
models the interaction of two biological species, one with local and one 
with nonlocal dispersal, that are competing for the same resource. We 
give a simple, concrete example of resources for which the equilibrium 
with only the local population becomes linearly unstable. In a sense, 
this example shows that nonlocal strategies can invade an environment 
in which purely local strategies are dominant at the beginning, provided 
that the resource is sufficiently sparse. 

Indeed, the example considered presents a high variance of the distribu¬ 
tion of the dispersal, thus suggesting that the shortage of resources and 
their unbalanced supply may be some of the basic environmental factors 
that favor nonlocal strategies. 
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1 Introduction 


The goal of this paper is to study the possible fate of a nonlocal diffusion strategy 
for a biological population in presence of a highly oscillating distribution of 
resource. 

The study of dispersal strategies and the comparison between local and 
nonlocal diffusive behaviors have recently attracted a great attention and several 
researches have been developed both in terms of experiments and from the 
purely mathematical point of view (see for instance [39l|29l[22l|34] and references 
therein). Remarkably, the phenomenon of possibly nonlocal hunting strategies 
has attracted also the attention of the mass-media, and related news can be 
found in popular newspapers and magazines (see e.g. m)- 

In this framework, even the distinction between local and nonlocal strategies 
is somehow a delicate issue and it is still not exactly clear in all situations what 
factors favor one behavior against the other. Of course, in general, as we know 
even from experience in our everyday life, it may be very difficult to deduce 
from overall principletfi] the optimal strategy to follow in each complex situation. 
Therefore, it is not surprising that the question of detecting the optimal strategy 
in a logistic mathematical model cannot have just a simple answer that is valid 
in every situation, and, concretely, very different dispersal strategies have been 
directly observed in nature. 

Detecting, analyzing and understanding the differences between diffusive 
strategies is therefore a difficult, but important, task in biology. One of the 
possible distinctions among the different strategies lies in rigorously defining 
the concept of “locality” (when a predator, roughly speaking, diffuses randomly 
in the neighborhood looking for an available prey) versus “nonlocality” (the 
short periods of hunting activity are followed by rather long journeys of the 
predator in the search for food). As expected, hunting strategies of predators 
are definitely influenced by the distribution of the resources. When the resources 
are “easily” available, it is conceivable that predators do not need to elaborate a 
nonlocal hunting strategy and indeed it can be more convenient not to drift too 
much to take advantage of the rather abundant resource in their neighborhood. 
Conversely, when the prey is sparse, it may be worth for predators to interchange 
the local hunting activity with suitable nonlocal travels in different possible 
regions. 

Of course, the more sophisticated the species involved in the hunt, the easier 
the latter phenomenon is expected to occur: namely, an intelligent species of 
preys will run away from the danger, thus making the distribution of resources 
for the predator sparse, and therefore making a nonlocal hunting strategy pos¬ 
sibly more favorable. However, in the model considered in this paper the re¬ 
source a is independent of the distribution of the populations, so this effect is 
not taken into consideration by the setting discussed here. 

^By “overall principles” we mean the availability of a general method, depending on the 
measurement of some parameters in the environment, which allows a population to choose an 
optimal strategy. We are referring to the impossibility of having a satisfactory and complete 
model for population dynamics, due to the complexity of the biological world. 
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It is also evident that the distinction between local or nonlocal strategy is 
a mathematical abstraction based on the consideration of different space/time 
scales: i.e., the ambient space that the population has at its disposal is not 
infinitely large in the real cases, and species cannot really perform discontinuous, 
nonlocal jumps. Nevertheless, a good mathematical model in which different 
scales are taken into account may furnish a justification for the diffusive strategy 
in a “large enough” environment in which the time scales of travel and hunting 
activities can be somehow distinguished in practice. 

We will try to give a rigorous mathematical framework to these naiVes consi¬ 
derations by showing the possible advantages of the long-jump dispersal strate¬ 
gies (i.e. the ones based on nonlocal diffusion) in regimes where the distribution 
of resources may be considerably different at different points of the ambient 
space. Not too surprisingly having in mind the concrete applications, we will 
use for this scope the mathematical framework of linearized systems and scaling 
properties of the eigenvalues, which take into account the stability property of 
equilibrium configurations. 

Our mathematical framework can be discussed as follows. Reaction-diffusion 
systems provide an effective continuous model for the biological problem of 
competition between different species. The typical example of local reaction- 
diffusion equation is 


ut = Am + {a — u)u in (0, T) x Q . (1-1) 

We study here the case of Dirichlet boundary conditions. Though other bound¬ 
ary conditions may be also taken into account to model different situations, 
our focus on the Dirichlet data is motivated by biological considerations (for 
instance, prescribing the solution to vanish outside a given domain corresponds 
to a confinement situation, for instance in a hostile environment). In this model, 
the environment is represented by the open bounded set 11 C M", with n ^ 2, 
and a heterogeneous resource cr : 11 —?► [0, -l-oo) is given (stationary in time). The 
growth of the population density u depends on a dispersal differential operator 
and on the reproductive rate of the population itself, which is proportional to 
the temporary availability of the resource {a — u). Dirichlet boundary conditions 
model a lethal environment for the population u outside the domain fl. 

A reaction-diffusion system involves at least two species, with distribution u 
and V, whose behavior is ruled by a reaction-diffusion equation like (11.111 . The 
two competing species differ for some special features: indeed, (HU) has to be 
modified in order to describe the foraging and reproductive habits of the species 
and further data concerning the environment. As it is customary in Adaptive 
Dynamics (see [19] and HD), the first step in the study of the evolution of a 
given feature is to single it out and then assume that the two populations differ 
for this feature only. For instance, in our mairH case the resident population has 
a local dispersal strategy and the mutant population has a nonlocal one. In [22] 

^Up to Section |2.3I we investigate the opposite situation, too, that is, when the resident 
population has a nonlocal dispersal strategy and the mutant population has a local one. 
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and m. one can find a comprehensive survey of the problem and of the standard 
approach in Adaptive Dynamics; many different features have been studied and 
compared in [3TJ [51 (TU] (different dispersal rates and genetic mutations), in [5T] 
(time-periodic sources) and in [TT1[T^|T3] (addition of a chemotactic component 
depending on the gradient of the resource). 

We are interested in the comparison of the dispersal strategies: in particular, 
we focus on the competition between a population with “standard” diffusion and 
a second population with nonlocal dispersal. Therefore, our model is 


Ut = Au + [a — {u + v)) u 

Vt = —{—A)^v + la — lu + v))v . 


( 1 . 2 ) 


At a discrete level, the “standard” assumption is that the motion of the po¬ 
pulation is governed by a random walk and this obviously leads to a Laplacian 
operator in the continuous model. Analogously, since our interest is focused on 
a second population with nonlocal dispersal, we adopt the fractional Laplacian 
operator as dispersal operator for the second distribution. The choice of such 
nonlocal diffusion operator is motivated by the fact that the fractional Laplacian 
has good stability properties in terms of the associated stochastic processes (it 
is the “continuous version” of the discrete motion governed by Levy flights, see 
e.g. [38] for a simple motivation and [3] for more advanced material), it possesses 
natural scaling features and it seems also to appear in real experiments (see 
e.g. [551115]). The present literature on the subject of nonlocal dispersal mostly 
considers convolution operators (see [HllHJIHl[T51|T5] ). In particular, in [55] . 
the model under investigation is 


Ut = fiAu -|- (cr — (m -I- v)) u 

vt = i' Jjy k (^) v{y) dy-v) + {a- {u + v)) v , 


(1.3) 


where p,, v are the dispersal rates of the two populations, respectively, and S is 
the dispersal distance of the second population. 

Of course, it is a delicate business to decide, in concrete situations, which 
models better describe the dispersion of a real biological population, and many 
nonlocal terms have been taken into account in order to comprise long-range 
effects. In general, we believe that fractional equations may be an important 
tool to further understand the complex problems arising in the mathematical 
modelization of biological species and we hope that the framework given in this 
paper can lead to a further development of the subject. 

In Section [5] we provide details and further explanations about the model 
considered here and some basic facts about the fractional Laplacian operator. 

We study the stability of a stationary solution (m, 0) of the aforementioned 
system, by means of a formal linearization at (w, 0), that we explain in Subsec- 
tion l2.31 The complete understanding of the global dynamics of a general system 
of diffusive and competing populations is beyond the scope of this paper and it 
seems, at first glance, very challenging from a mathematical point of view, since 
a variety of possible situations may occur. Nevertheless, let us stress that even 


4 


the analysis of the stability of a stationary solution (also called “invasibility” 
analysis) is interesting and meaningful from an evolutionary point of view, as it 
is suggested in the principal literature in Adaptive Dynamics (again, see |19|b 
In fact, a small perturbation around (u, 0) mirrors the occurrence of a genetic 
mutation in the first population, involving the dispersal strategy. At (u, 0) the 
hrst population benefits from an equilibrium state, while the second one does 
not even exist. Then a small portion of the first population (with density u) 
undergoes a genetic mutation, which starts a second population (with very small 
density v) which competes for the resource with the former. Of course, the ge¬ 
netic mutation of this theoretical experiment involves only the hunting/dispersal 
strategy, passing from a local to a nonlocal one. In this context, the expected 
outcome of the analysis of the stationary solution is, in most of the cases ex¬ 
perienced in practice, stability, that is, the second population does not find the 
right conditions to evolve and it gets rapidly extinguished. On the contrary, 
(even partial) instability of these type of equilibria is rather surprising and in¬ 
teresting, since in this case the new dispersal strategy is convenient enough to 
allow a short term survival of the second species and to provide a situation of 
coexistence of two different populations. 

The core of this paper is Section O where we show how the stability of (u,0) 
(namely, the sign of the eigenvalues associated with the linearized system) de¬ 
pends on the distribution of the resource a. In particular, we will show that if 
a certain relationship between the variation of a and the fractional Poincare- 
Sobolev constant in fi is fulfilled (see Definition 13.1|) . then the linearized system 
has a positive eigenvalue and (u, 0) is unstable. It is transparent from Dehni- 
tion O that the distributions leading to instability of (u, 0) (and suggesting 
convenience of a nonlocal dispersal strategy) are those with a “huge variation”. 
The last part of Section |3] is devoted to show that such a distribution a may 
occur. Summarizing, the result that states that the local dispersive strategy 
may become unstable in presence of a new population endowed with nonlocal 
diffusive strategies can be formally stated as follows: 


Theorem 1.1. Let fi C K" be an open subset o/R" with Lipschitz boundary 
and let s £ (0,1). There exist bounded functions cr : D —>■ [0,-|-oo) and u : 
ft —> [0, -boo) such that (u, v) := (u, 0) is a linearly unstable equilibrium for the 
system 


{ Ut = Am -|-(<t — (m -b m))m in ft 
Vt = —(— A)®v -b(cr — (m -b v))v in ft 
u = 0 on dft 

V = 0 m M" \ D . 

More precisely, the function u is a solution of 

f Au(x) + (a'(x) — u{x))u{x) =0 in ft, 

1 M = 0 on dft, 


(1.4) 


(1.5) 


and the linearization of system in m at (m, 0) has a negative and a positive 
eigenvalue. 
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The existence of distributions of resource a which support the phenomenon 
described in Theorem 11.11 is motivated by concrete models. To be precise, in 
this paper we examine two particular cases: 

1. a rescaled resource u\ on a sufficiently small domain fix in Section [3. II 

2. a large multiple of the characteristic function of a ball in Section 

We remark that Theorem 11.11 states that u is a linearly stable solution of 
the autonomous, scalar Fisher-KPP equation in (II.5L but (u,0) is a linearly 
unstable equilibrium for the system in (11.41) . More explicitly, the positive eigen¬ 
value of the linearized system takes into account the fact that if the density 
of the first population undergoes a small variation without the appearance of 
the second species, then the system has the tendency to return to the original 
position. Conversely, the negative eigenvalue shows that if a second population 
appears, then the system does not go back to the original situation, and the 
second species has indeed chances to survive and colonize the environment. 

It is interesting to contrast this result with those obtained in [35] : when the 
dispersal rates are equajl, that is fj, = v, if the dispersal distance (5 in (II. 3p is 
sufficiently small, then v can invade the local population u but u cannot invade 
the nonlocal population v. This suggests, as a general principle, that the smaller 
spreader may be favored by evolution, especially in hostile environments. Our 
approach is rather different: in some sense, we consider the dispersal distance 
as already fixed (in the definition of fractional Laplacian) and we investigate 
the dependence of the possibility of an invasion based on the availability of the 
resource a. Our attempt is to put in evidence the role of the environment (more 
precisely, the fact that Definition 13.II is fulfilled) in the selection of the dispersal 
strategy. 

Roughly speaking, the condition (in Definition 13. II) which allows the instabi¬ 
lity of the system records the fact that the first population, with local diffusion, 
cannot saturate the given resource and leaves enough “leftovers” for the second 
species to survive. 

In this sense, a natural question is to determine whether a population ex¬ 
hausts the resource. For this, as a second result, we provide an example of a 
purely nonlocal phenomenon in population modeling. We show that, fixed any 
arbitrarily small e > 0 and given any resource a G C^{Bi, [0, -boo)), there exists 
a resource G C^(Ri, [0,-boo)) that is e-close to cr in the norm of 

a radius > 1 and a function which vanishes outside which is s- 

harmonic in Bi, which equals to in Bi and which therefore satisfies 

(-A)®Me = (CTe - Ue)Ue in Bi. 

That is, up to an arbitrarily small error, a nonlocal population can locally adapt 
to any given resource (provided that the density of the population is artificially 
and appropriately regulated in a suitable region). The formal statement of this 
result goes as follows. 

®In our model, we do not even take into account different dispersal rates fi and u. 
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Theorem 1.2. Let fc € N and a G C^{Bi, [0, +c»)). Fix e > 0. Then there 
exists CTg G with 

Ik-o'ellcnsi) ^ ^ (1-6) 

and there exist > 1 and G C^{Bi) fl C"*(R") such that 

Ue{x) = <7e{x) for any x G Bi 
{—Ayui;{x) = 0 for any x G Bi 
Ue{x) = 0 for any x G M" \ B^^^ 

In particular 

{-AYue{x) = {ae{x) - Ue{x))ue{x) for any X G Bi . (1-10) 

It is worth mentioning that Theorem 1 1.2 1 heavily relies on the nonlocal fea¬ 
ture of the equation and it does not have any local counterpart (this will be 
clearly explained in Section SI). Let us stress the fact that Theorem Q] does 
not prove (and cannot prove, since this would be false in general) that a non¬ 
local population always exhausts completely the resource, since a small error e 
has to be taken into account. In a sense, the solution given by Theorem 11.21 
is different than the original one, since it does not attain the homogeneous 
Dirichlet boundary datum: it has prescribed, non-homogeneous (but compactly 
supported) Dirichlet boundary datum outside the strategic region in which the 
equation is satisfiec0. 

We observe that Theorem O has important (though socially not embrace- 
able!) practical consequences. For instance, a given population may have a 
strong intention to consume all the given resource in a region of particular 
strategic importance (say, a region contained in the ball Bi in our example). 
Indeed, in concrete cases, this strategic area might be favorable for generating 
a new competing species, or might be easily accessible by a similar population 
coming from abroad which can be considered dangerous or undesired by the 
local population, and the possible leftover of the resource might obviously fa¬ 
vor the newcomers. For these reasons, a “socially conservative” (and rather 
unkind!) population may wish to avoid to leave available resources in strategic 
regions which can be used by unwanted competitors. 

The result in Theorem 11.21 savs. roughly speaking, that in this case, a non¬ 
local population is able to find a suitable, somehow “artificial”, distribution of 
population far away in order to consume the resource in the strategic region and 
thus penalize the newcomers (viceversa, a local population cannot do that). 

Notice that this suitable distribution of the conservative population may 
require a modification of the conditions far away: indeed the supporting ball 
Bfi^ ^ in Theorem 11.21 may become larger and larger for small e: that is, in a 
sense, the conservative population may need to change its plan “close to infinity” 

^The fact that the Dirichlet boundary condition is not homogeneous reflects mathematically 
the practical condition of performing an effective distribution plan for the population outside 
the strategic region. 


(1.7) 

( 1 . 8 ) 
(1.9) 
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in order to consume more efficiently the inner resource (in this sense, a “global 
plan” for the population distribution is in order, and it is indeed conceivable 
that an optimal use of resources may involve strategic plans on the distribution 
of the population in the large). 

The rest of this paper is organized as follows. In Section[2]we recall the basic 
notation about the population dynamics model that we study. The linearized 
dynamics of the system is then analyzed in Section [3 where we will also give 
two examples that establish Theorem 11.11 Finally, in Section U we will prove 
Theorem ll.2l and show that it is a new phenomenon, which only arises in nonlocal 
dispersion models. 


2 Biological models and mathematical tools 

2.1 Population dynamics 

Let us denote hy u,v : [0,T) x fl — [0,+oo) the densities of two species coex¬ 
isting in the same domain fl and competing for a common resource tr : —>■ R. 

Here and in the rest of the paper we consider as a domain an open, bounded set 
n C R" with Lipschitz boundary 9H. The resource a belongs to the space of 
measurable, essentially bounded functions We study the linear stability 

of a stationary solution of the reaction-diffusion system with Dirichlet boundary 
conditions 


{ Ut = Art -I- (cr — (it -|- u)) u 

vt = -{-AYv -I- (cr - (it -I- v)) V 
u{t,-) = 0 
= 0 


in [o,r) X n 
in [o,r) X n 

on 912, Vt S [0,T) 
in R’"\H, VtG [0,r). 


( 2 . 1 ) 


For this, we perform a formal linearization around a stationary point (it, 0) of 
(ED and then we focus only on the corresponding linearized system, that is 

{ — All = (cr — 2u)u — uv in H 

(—A)®i; = (cr — u)v in H 

It = 0 on dQ. 

i; = 0 in R” \ H . 

Remark 2.1. Though the global dynamics is beyond the scope of this pa¬ 
per, we recall that there is a detailed and specialized literature about the well- 
posedness of the initial value problem associated to (12.11) with it(0,x) = uo{x) 
and ^(Ojx) = vo(x) for some given functions uo,vo (see [35] for the well- 
posedness of the problem in for instance). The analysis of the global 

dynamics of semilinear parabolic systems is performed through the theory of 
Monotone Dynamical Systems (see [26], [2^ and [36]). For the fractional Fisher- 
KPP equation one should see [5] and [1], for instance. See also Section O for 
further remarks. 


Before focusing on the aforementioned linearized system, let us recall some 
useful definitions and facts about the pseudodifferential operator (—A)® that is 
involved in (12.11) . 


2.2 The nonlocal dispersive strategy and the fractional 
Laplacian 

Consider an open set C M" and s G (0,1), the Gagliardo seininorm of a 
measurable function u is defined as 




/ / 

JR^ JR^' 


\u{x) - 
|a; - ?/|"+2s 


dx dy 


The fractional Sobolev space that we denote here is the linear set con¬ 

taining all the measurable functions u : R" —?► R such that: 

• I|m||l 2 (q) < -hoo, 

• < - 1 - 00 , and 

• u{x) = 0 for a.e. x G R" \ fl. 

The Gagliardo seminorm is naturally related to the fractional Laplacian, since 


(-A)®m(x) := 


T{n/2 + s) \ f u{x) - u{y) 

^2s+n/2Y(^_,)) |x - y|"+2® 


dy, 


where T is the Euler’s function. For an introduction to the fractional Laplacian 
and the fractional Sobolev spaces see for instance [TH]. In our framework, the 
scalar version of (EH), that is 


vt = —(—A)®u -I- (cr — v)v , 


is known as Fisher-KPP equation with fractional diffusion and for the many 
established results one can see, for instance, [5] and m- 

In this section we summarize the results needed in this paper only. 

Theorem 2.2 (Fractional Poincare-Sobolev embedding theorem). Fix s G (0,1) 
and an open bounded set O C R" with Lipschitz boundary. There exists a positive 
constant Cj = (7)1(3,0) such that 

y (j) G Hq{TI) , ||^||i2(Q) ^ (7)j[(/)]^s(Rn) . (2.2) 

This means that HQ{n) is continuously embedded in 

Proof. We give the proof, which is of classical flavor, for the facility of the reader. 
We argue by contradiction, supposing that there exists a sequence (j)k G Hgi^) 
such that ||(()fc||L 2 (Q) ^ We define 


WMl^q) 
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Then G and 



|' 0 fc(a^) - - 0 fc(y)P 

|a; - ?/|"+ 2 '> 





\\ 4 >k\\L'^(Q) ^ k 


Therefore 



Itl^kjx) - 'ipkiy)\'^ 
|a; - ?/|"+ 2 * 




Also, ||V^fc||L 2 (n) = 1 . Therefore, by compactness (see e.g. Theorem 7.1 in [T 8 ] . 
used here with p = q = 2), we obtain that, up to a subsequence, ipk converges 
to some tp in L^{Vl) and a.e. in 17. Defining '4>{x) := 0 for any x G K” \ 17, we 
have that ■0^ = '0 = 0 a.e. in R” \ 17, and consequently ipk converges to ip a.e. 
in K’". 

Thus, by taking the limit in (12.31) and using Fatou’s Lemma, 



\ip{x) - ip{y)\^ 
|a; - ?/|"+ 2 s 


dx dy ^ lim inf 

k—^-\-oc 


^ lim inf 7 ^ = 0 . 
k — K 


| 0 fc(a:) - - 0 fc(l/)P 

\x - y|’^+ 2 s 


dx dy 


Accordingly ip must be constant in R" and therefore identically equal to zero 
(up to sets of null measure). This implies that 


1 = lim || 0 fe||i 2 (o) = lim || 0 fe - 0 ||L 2 (n) = 0 . 

k—¥-\-oo k—^-\-(X) 


This is a contradiction and it proves the desired result. □ 

In the following, we will always assume (^[((s,!!) to be the sharp constant 
such that ( 1 ^ holds, namely 


C7\s,n) = ^ inf 








2 

L 2 (n) 


= inf 




4ieHg(n) ||0||^2mn 
<^.^0 ^ 


(2.4) 


Remark 2.3. If r > 0 and 0 G Hq{Bi), one can consider the rescaled func¬ 
tion (pr{x) := r“"/^0(a;/r). Then 0^ vanishes a.e. outside Br- Moreover, 
ll0r|U2(Rn) = ||0||i2(Rn) uud [0r]ff-(R") = r-'’[0]j^.(Rn). Accordingly, 

Q(s,B,) =r 2 ^Q(s,Ri). 


2.3 Linearization of the system 

Let 17 C R" and a G L°°{il) be as in Section [51 Our purpose is a qualitative 
study of an equilibrium state of the following system 

( ut= Au +{(T-{u + v))u , , 

( Vt =—{—AYv+ {a— {u + v))v ^ ' ' 
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More precisely, we look for an equilibrium state of the form {u, 0) with u G 
i?g (n) and u > 0. 

For the sake of completeness, we also investigate the existence of an equilib¬ 
rium state of the form (0,1;), with v G v ^ 0. The linearization of 

(|2.5D at (0,f;) and further conclusions are postponed to Section SI 

Definition 2.4. Given a bounded function cr : O —>■ [0-1- oo), we say that a 
satisfies a reverse Poincare-Sobolev condition if 

sup / a{x)u{x)^ dx — / \Vu\^dx>0. (2.6) 

u&h^{Q)Jq Jq 

Furthermore, cr satisfies a reverse fractional Poincare-Sobolev condition with 
parameter s if 

sup / a{x)v{x)'^ dx - [vffjsmn) > 0 . (2.7) 

Jn 

In order to make computations easier, we give a sufficient condition that 
ensures (EH). 

Lemma 2.5. Let Ai (f2) be the first eigenvalue of the Laplacian in Q with Dirich- 
let boundary condition and let G Hq{D) be the corresponding eigenfunction. 

If 

Ai(n) j (j)i{x)^ dx < I a{x)cj)i{x)^ dx , (2.8) 

JQ Jn 

then the reverse Poincare-Sobolev condition in (1^ is satisfied. 

Proof. By construction 

f —A(/)i = Ai(r2)(/)i in fl 
1 (()i = 0 on dLl , 


and so, by (12.81) . 

j |V()>i(x)|^ dx = Ai(0) j (j)i{x)'^ dx < j (j{x)(j)i{xY dx 

Jq Jq Jq, 

which proves (1^ . □ 

Remark 2.6. It is worth noticing that condition (12.8|) is satisfied, for a fixed 
domain LI, for any resource a that is sufficiently large in an open subset of LI. 
Hence, fixed LI, there are many examples of smooth resources satisfying (12.8|) 
and therefore (12.61) . 

Remark 2.7. We also observe that the converse of Lemma [2.51 does not hold 
true, i.e. the reverse Poincare-Sobolev condition in EH does not necessarily 
imply (12.8|) : as an example, one may consider LI = (0, tt), a{x) = e“^®'^^°X(o,e)(2;) 
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and u{x) = with £ > 0 suitably small. Then u € and (12.61) holds 

true, since 


/ a(x)u(x)^ dx — f |Vupdx = £ f dx — ^ [ x dx 

Jn Jn Jo 9 Jo 


3 

“ 7^ 


- 29 / 10 ^ 7/3 _ 477 ^/^ ^ 


On the other hand, in this case 4>i{x) = sina:, Ai(r2) = 1, and 
Ai(r2) / <l)i{x)‘^ dx — / a{x)(l)i{x)'^ dx 

Jq, Jq 

= j sin^ xdx — £-29/10 f g-j^2 ^ 

Jo Jo 

TT £-29/10 ^ 

=-(£ — Sin£COS£) . 

2 2 ^ ’ 

Thus, since, by a Taylor expansion, 

sin£cos£ = (£ + 0(£^))(1 + 0{e^)) = e + 0{e^) 


it follows that 

£-29/10^£ _ sin£COS£) = 0(£^/^°) 

and so 

Ai(0) f (f)i{xY dx — f a{x)cj)i{x)^ dx = ^ > 0, 

Jn Jn 2 

which shows that (j)i does not satisfy (12.8|) . 

The fractional equivalent of Lemma [23] is stated in the following lemma. We 
omit its proof, which would be a repetition of the proof of Lemma 12.51 

Lemma 2.8. Let (Js G Hoi^) ® minimizer for the Rayleigh quotient (|2.4I) . If 


Cj ^(s,f7)||(^s|||2(n) < J^<j{x)(j)s{xf dx , 

then the reverse fractional Poincare-Sobolev condition with parameter s in (1371) 
is satisfied. 

Remark 2.9. As we noticed in Remark 12.61 for the local analog of the reverse 
fractional Poincare-Sobolev condition, if the resource a is sufficiently abundant 
in the domain LI, then condition (12.71) is satisfied. 

The reverse Poincare-Sobolev condition in (1731) is a useful tool to obtain 
non-trivial solution of the local stationary equation, as stated in the following 
result. 
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Theorem 2.10. Consider s € (0,1] and a bounded function tr : 12 —>■ [0,+c») 
satisfying either 

(1) the reverse Poincare-Sobolev condition in (12.61) (when s = 1) or 

(2) the reverse fractional Poincare-Sobolev condition in (12.71) (when s < 1/ 

Then there exists a non-trivial, non-negative function u € Hq{Q) (i.e. u ^ 0 
and u^O) satisfying 


f (—A)^u(x) = (cr(x) — u(x))u(x) 
[ u = 0 


in 17 
on dQ. 


(2.9) 


Proof. The proof is a minimization argument, based on coercivity and energy 
methods. Though the idea of the proof is rather standard, see e.g. [S], we 
provide the necessary details for the facility of the reader. The proof also gives 
us the possibility of a comparison between local and nonlocal case. 

First, we prove the theorem for the local case s = 1 and then we provide the 
suitable changes in order to prove the nonlocal case, too. 

(1) Consider the following energy 


E{u) 




3 


defined on i7Q(17). Notice that the Euler-Lagrange equation for E gives 


—Au = {a — |u|)m . 

We show that the energy E is coercive in i7Q(17), that is 

E{u) —>■ +00 as ll'ull—>■ +00 . (2-10) 

For this, we use the Young inequality with exponents 3/2 and 3 to see 
that, for any a, 6 ^ 0, 


, 2 3 1,., 

ab ^ -a^ + -6^ . 

In particular, taking a := and b := 

that 


a 


II II 


2-^/^lklU °°(a)i we obtain 

lklli»(o) 

6 


hence 


+ — ^ -Co, 


for some cq > 0 independent of u. Accordingly, 


JO. ^ 
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that establishes (12.101) . 


As a consequence of (I2.10L we have that E has a global minimum u G 
Hglfl), satisfying 

—Au = (a — \u\)u . 

Since u is a minimum, then u := |?l| is a minimum too, because E(u) = 
E{\u\). Thus we can consider a non-negative function u ^ 0 satisfying 

—Am = (a — u)u . 


We conclude the proof by showing that condition (12.6|) guarantees that 
E{u) < 0 and then u ^ 0. By (12.61) . there exists a function u G 
with 


/ a{x)u{x)^ dx — / \Vu\^dx>0. 

JQ Jn 


By density, we can suppose that u G (7^(12). For every e > 0 we can 
rewrite the energy E evaluated at eu as 


E{eu) = 



hence E{u) < E(eu) < 0 provided e: is small enough. 
(2) The energy 


Es{v) ■■= 


R"xR" 


\v{x) - v{y)\'^ 
|a; — ?/|"+2s 






is well definecd and coercive in iJg (O) and the proof is the same as in the 
local case. Moreover, the Euler-Lagrange equation for Eg is 

(—A)®?; = (ct — |m|)u . 

Consequently, Eg has a global minimum v and v = |U| is a minimum, too, 
because 


Eg{y) = 


|M(a;)| - |r;(y)|p 






R^xR” |x — yp+2® 
v[x) -v{y)\'^ 
\x - 


I ”—+ 3 






dxdy- / (j—+ '—= Eg{v). 
Jo. ^ o 


As we proved in part (1), condition (12.71) ensures that Eg{v) < 0 and thus 

M ^ 0. □ 


^For the sake of simplicity, we omit the multiplicative normalization constants. 
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The result in Theorem l2.10l and several variations of it are rather of classical 


flavor: with slightly different assumptions on a (take, for instance, ct > 0 in 
n) and a branching condition matching (12.61) for the existence of non-trivial 
solutions, it can be found in [5] and in [3]. 

In view of Theorem 12.101 and comparing with Remark 12.61 and Remark 12.91 
we obtain the fact that the richer the environment is, the easier the survival of 
a population. This fact, which matches the intuition, finds a detailed quantifi¬ 
cation in the following observation. 

Remark 2.11. Since the reverse (fractional) Poincare-Sobolev inequalities (12.61) 
and (12.71) seem to play a symmetric role in Theorem 12.101 let us compare them 
more carefully. In some sense, the Dirichlet boundary conditions being equal 
for (12.91) when s = 1 and when s < 1, the nonlocal population has an advantage 
when the diameter of the domain tends to 0. 

More precisely, as we remarked in Remark 12.61 a resource a needs to be 
sufficiently large in order to meet (12.81) . which implies (12.61) . How large should 
a be is proportional to the first eigenvalue of the Laplacian Ai(H) = C^{1, 17)“^. 
Now, if H = Hr, we observe that 


1 



as r —0 , 


that is, the environment becomes more and more lethal for the local population, 
because (12.81) is very difficult to satisfy. The situation is milder for the nonlocal 
population, because, thanks to Remark [3131 



as r 0 . 


This means that the criticality of the domain size is slower to prevail on a 
nonlocal population. 

In this sense, the lethal property of the boundary (as described by the ho¬ 
mogeneous Dirichlet datum outside the domain) has a different influence on 
local and nonlocal populations, depending on the scale of the domain. For in¬ 
stance, for small balls (when 1/r ^ l/f®), nonlocal populations are favored. 
Conversely, for large balls (when 1/r <C l/r®), local populations are favored 
(heuristically, because the local diffusion has little chance to reach the deadly 
boundary). 

In the perspective of an applied analysis, one can find explicit fractional 
Sobolev constants in m- 

In the remaining part of this section, we focus on the local case, that is, on 
the properties of a stationary point for the system (1331) of type (11,0). This is 
motivated by the evolutionary point of view of studying the effect of the advent 
of a new population (this approach is indeed often adopted in the literature, see 
e.g. [21], [24]). Of course, we think that it would be also an interesting problem 
to investigate the cases of a dominant nonlocal population (corresponding to a 
stationary point of type {0,v)) and of the possible coexistence of two different 
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populations, namely provide concrete assumptions on the resources and the 
domains that allow the existence of equilibria {u*,v*) with both u* and v* 
nontrivial. 

Of course, an easier approach to the existence of mixed states {u*,v*) may 
be taken by studying the case of different resources in the two logistic equations, 
but we do not address this problem in the present paper. 

Remark 2.12. As a byproduct of the proof of Theorem l2.101 we have that the 
solution found is an energy minimizer. That is, if u is the solution obtained in 
Theorem 12.101 then E{u + eu) ^ E{u), for any u G Accordingly, the 

map 

e !->■ £{e) := E{u + eu) 
attains its minimum at e: = 0 and therefore 

0^£:"(0)= / \Vu\^ - au^+ 2uu^ dx. (2.11) 

Jn 

In particular, the solution is linearly stable, i.e. the second derivative of the 
energy is a positive quadratic form. 

The energy functional is quite useful to capture the stability of the pure 
states, such as the ones of the type (u,0). For related approaches, also based 
on the linearization of semilinear systems, see e.g. [5]. Also, it is useful to recall 
that the population u cannot beat the resource cr, as stated in the following 
result: 

Lemma 2.13. Consider a bounded function cr : fl —> [0,+oo) and a non¬ 
negative solution u € of (12.91) . Then u{x) ^ ||cr||ioo('Q), for any x G £1. 

Proof. Let 0 := ||cr||2,oo(n). We test equation (12.91) against v := max{M — 0, 0} 
and we see that 

/ |Vup = / Vit • Vu = / {a — u)uv = / (a — u)u(u — Q). 

Jq Jq Jq J {u^Q} 

Now observe that, in {{t ^ 0}, we have a — u^Q — u^O, which shows that 


|Vwp 0. 


Accordingly, v vanishes identically and so m ^ 0. □ 

Corollary 2.14. Consider a bounded function cr : 17 —>• [0, +oo) and a non¬ 
negative solution u G iLg(f7) of (12.9|) . Then u is continuous inside 17. 

Proof. One defines 0 := ||cr||ioo(f2) and tests equation (12.91) against v := max{u— 
0, 0} to obtain the desired result (see e.g. [8]). □ 
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From now on, we focus on the stability of the system around the stationary 
point (u, 0), where the distribution of resources a satisfies (12.61) and u G Hq(Q) 
is a non-trivial, non-negative solution of (12.91) . 

The linearization of the system (12.51) at (u, 0) gives, as a result, the linear 
operator 


L(u,o){u,v) 


A -I- (cr — 2'ii) —u \ / u 

0 -(-A)®-I-(cr - it) y V ^ 

Am -|- (cr — 2m)m — uv \ 

— (— A)®m -I-(cr — m)m J ’ 


( 2 . 12 ) 


for any (m,m) G i(io(fl) x Fig (12). The associated quadratic form, with respect 
to the duality in Fig (12) x ii/g)!!), is 


/ {cr-2u)u^-uuv+{a-u)v'^ dx, (2.13) 
J Q 

for any (m, v) G iifg (12) x iJQ(12). From the triangular form of T(s,o)i the relevant 
information is concentrated on the signs of the principal eigenvalues of the pseu¬ 
dodifferential operators on the diagonal of (j2.12|) . In this spirit, we first point 
out that the direction (u, 0) is always linearly stable. This is pretty obvious if 
we think at the biological model, since (m, 0) is the stationary configuration of 
just one population, and slightly and proportionally modifying the density of 
this population without letting any new population come into the environment 
should not drive the system too far from the previous equilibrium. The formal 
statement goes as follows: 

Lemma 2.15. As long as there exists a solution u for (EH), we have that 

<3(fi.o)(u:0) < 0- 

Proof. By testing (12.91) against m, we obtain that 


As a consequence. 


Jffi(R") 



m)^m^ dx. 


Q(ii,o)(M:0) =-[w]ffi(R") + / {a - 2u)v? dx = - / vA dx. 

Jn Jn 


The latter term is strictly negative, thanks to Theorem 12.101 and so we obtain 
the desired result. □ 


We point out that Lemma[2T5]is a particular case of a more general stability 
result. Namely, the stationary configuration (m, 0), which corresponds to the 
local population colonizing the whole of the environment, is also linearly stable 
with respect to all the perturbations in which only the the density of the local 
species varies (i.e. the possible source of instability in this setting may only 
come from the advent of a nonlocal population). The formal result goes as 
follows: 
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Lemma 2.16. As long as there exists a solution u for (EH), we have that 


Q{u,0) 0 ) ^ 0 

for any u £ 

This lemma is well-known, due to the variational characterization of the 
associated eigenvalue problem. We include a proof for the convenience of the 
reader. 

Proof. From (I2.13|) . 

Q{u,o){.'U',0) = — / \Vu\'^dx+ / {a — 2u)u^dx, 

JQ. JQ 

hence the claim follows from (12.111) . □ 

In view of Lemma 12.151 we obtain that a good way to detect the possible 
linear instability of the point (u, 0) is to rely upon the perturbations of the 
form (0, v), i.e. in the possible advent of a new population with different diffusive 
strategy. The purpose of the next section is therefore to understand when it is 
possible to obtain that 

Q(u,0) p 

for a suitable choice of u* G iJg(O). 


3 Linear instability 


Our aim in this section is to enlighten the connection between the distribution 
of resources a and the possible instability of the system, which would suggest 
some convenience in a nonlocal dispersal strategy of the second species v. For 
this, we introduce the following notation: 

Definition 3.1. Let a : —>■ [0,-|-c») satisfy the reverse Poincare-Sobolev 

condition of Definition 12.41 Let u ^ 0 be a non-trivial solution of the non¬ 
linear equation (12.9L provided by Theorem l2.10l We say that the pair (ct, ft) is 
mismatched in D if there exists ccg G D and r > 0 with Br(xo) C fl and 


inf ((T(a:) - u(a:)) > 

xGBr(xo) 


1 

Cjj ( 5 , Br (^ 0 ) ) 


(3.1) 


In this formula, the constant Cjj(s, Br(xo)) is the sharp fractional Poincare- 
Sobolev constant with respect to the ball Br.{xo) provided by Theorem 12.21 

Roughly speaking, condition EH says that the solution u is not capable to 
exhaust the whole of the resource a in the whole of the domain: that is, at least, 
in the region Br{xo), the population does not manage to take advantage of all 
the resource at its disposal and there is at least a quantity Cj|(s, as a 

leftover. 
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In Subsection 3.2 we will see an example of mismatching (cr, ft) and it will 
be clear in that case that the mismatch condition depends basically on a only. 

In our setting, condition (13.11) is sufficient to ensure linear instability, as 
given by the following result. 

Proposition 3.2. // the mismatch condition in en) is satisfied, then there 
exists Vi, € such that Q(iip)(0,i;*) > 0. 

Proof. By (12.41) and en), we know that there exists Xq G and r > 0 such 
that 

Br{xo) C fl (3.2) 

and 


inf ((j{x) - u{x)) > ^^ 

xeB,.{xo)^ ^ C(,{s,Br{xo)) 


inf 

<j>eH‘iBr{xo)) 

^^0 


[4‘]'h‘( m.n') 
W'^Wl^iBr-ixo)) 


As a consequence, there exists z;* £ HQ{Br{xf)) such that z;* ^ 0 and 


• r II \ 

mf [a(x) - u[x)) > J —^2-. 

xeBAxo)^ \\v4hiBAxo)) 


(3.3) 


Now notice that ||z;*||i2(s^(3,jj)) = ||z;*||i2(f2) and z;* vanishes a.e. outside fl, 
thanks to (13.21) . This gives that z;* G Moreover, by (12.131) and (13.31) . 


Q{u,o)iO,Vi,) = -[z;*]^. 


' Br(xo) 


(a — u)v^ dx 


V 


> — z;. 


*\H‘‘ 


+ 








L^{Bi{xo)) JBr(xo) 


' dx = 0. 


which gives the desired result. □ 

Remark 3.3. Proposition 13.21 proves the linear instability of the point (zz, 0) 
with respect to perturbation of the type (0,z;*) (compare with the theory of 
Monotone Dynamical Systems in [36] or see El)- 
Indeed, 

Q{u,o){0,Vi,) = ||i;*|||2A(D) , 

where A(D) is the principal eigenvalue of the linear pseudodifferential operator 
— (—A)®+((t— zz) (see the characterization of the principal eigenvalue by Rayleigh 
quotient in |4]), thus T(i,o) has a negative eigenvalue and a positive one and the 
stability of a stationary state is determined by the spectrum of the linearization 
(for this general principle see [55]b 

Heuristically, this can be understood as follows: by formally plugging (zz, v) = 
(zz, 0) + £(0,z;*) + o(e) into (11.41) we obtain 


Vt = —(—A)®z; + (cr — (zz + v))v = —£(—A)®z;* + e{a — zz — £z;*)z;* + o{e). 
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Thus, since Vt = edtVi, + o(e), we formally obtain 


dtVi, = -(-A)'*'u* + (o- - u)u* + o(l). 


Hence 

5ilk*lli2(Rn) = 2 / u*atu*da; = (3 (£,,o)(0,'c*)+ o(l), 

jR" 

which is positive by Proposition [321 

Therefore.ProDOsition 1321 states that the size of the new population (mea¬ 
sured in the T^-norm) has chances to increase (at least for short times). 

These type of linearization arguments in the neighborhood of equilibria that 
correspond to only one biological species are widely used in Adaptive Dynamics, 
see for instance [13], [33], 130] and the references therein. 

The rest of this section is devoted to show that the assumptions of Proposi¬ 
tion 13.21 hold for some cr : D —>■ R. 

3.1 Rescaling arguments 

We propose here a rather simple rescaling argument which gives the existence 
of a domain and a distribution of resources <t\ satisfying the assumptions 
in Proposition 13.21 The main drawback of this argument is the fact that the 
domain changes with the parameter. On the other side, it is immediately 
evident that the resource a\ leads to instability at (uajO) when it starts being 
sparse and far from being homogeneous. 

We consider here a smooth function cr : D —>■ [0, -l-oo) satisfying the reverse 
Poincare-Sobolev condition in (|2.6I1 (recall Remark |2.61) and the corresponding 
stationary solution u given by Theorem 12.101 We see that, in this case, the 
population u does not exhaust the resource a in the whole of D. More precisely, 
we have: 

Lemma 3.4. Let cr : D — >■ [0, -|-oo) be a smooth function satisfying the reverse 
Poincare-Sobolev condition in (1^ and let u be the corresponding stationary 
solution given by Theorem \2.10[ 

Then there exist xq S D, r > 0 and cq > 0 such that Br{xo) C D and 

o'{x) — u{x) ^ Co 


for any x € Br{xo). 

Proof. By testing (12.91) against u, we obtain that 

0< / \Wu{x)\'^ dx = / {a{x) — u[x))v?{x) dx. 

Jq Jq 

This implies that there exists xq € such that cr(xo) — u(xo) > 0. The desired 
result follows from the continuity of u given by Corollarv l2.14l □ 
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In the notation of Lemma 13.41 by possibly translating the domain, we can 
assume that xq = 0, and so 


a — u ^ cq > 0 in Br- 


(3.4) 


Then we consider the family of rescaled domains 

ilx ■■= {X~^y : 2/ e fl} 

and rescaled functions 

a\(x) := Xa{'/Xx), Wx G fix 


with X ^ 1. Then 

ux(x) := Xu{^/Xx ), yX G ilx 

is a positive stationary solution for the equation (12.9|) with resource cta, since 


{Aux + {ax — ux)ux) (x) = {X^Au + X^(a — u)u) {VXx) = 0 , Vx G fix ■ 

Proposition 3.5. There exists A ^ 1 such that, for every A ^ A, the pair 
{ax,ux) is mismatched in the corresponding domain Qx, according to Defini- 
tion \S.l\ 

Proof. We take rx := X~^r. By (13.4|) . 


inf (tTA(a;) - MA(a;)) 
|x|<rx 


inf ^ A {a{'/Xx) — ux{'/Xx)) 

|a;|<A~ S ’r 


= inf A (cr(y) - ux{y)) ^ cqA. 

|y|<r 


On the other hand, by Remark 12.31 

Ci{s,Br,) = rl^Ciis,B,) = X-^r^^Ci{s,B,). 


(3.5) 


By comparing this with (13.5L we conclude that 


Ci{s,Br^{xo))’ 


provided that 


A> (cor2«Q(s,Bi)) 


□ 


From Propositions 13.51 and 13.21 we obtain that there exists G Bg(fix) 
such that (5(iLj,,o)(0 ;'*^*.a) > 0, as long as A is large enough, hence (ua, 0) is 
linearly unstable. 

This is a first example that shows the validity of Theorem 11.11 (a different 
one will be constructed in the remaining part of this paper). It is worth pointing 
out that the condition that A is large translates into the fact that the domain Qx 
is small and the resource ax is very unevenly distributed. In some sense, the 
nonlocal diffusion may allow the population to take advantage of the small region 
in which the resource is abundant, while a less diffusive population may starve 
in the portion of the environment with limited resource. 
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3.2 Branching arguments 

In this subsection we focus on a particular family of distributions, indeed we 
assume C and 


<yr{x) := TXBr{,xo){x) 


T X G Br{xo) 
0 X ^ Br{xo) 


We show that there exist r, r > 0 such that the assumptions of Proposition 
[SJ hold. First of all we have to deal with with Definition [211 which located a 
branching point for solutions of (12.91) . For this, for any r G K, xq £ K”, r > 0, 
such that Br{xo) C D, we introduce the quantity 


e{T,xo,r) := sup t — / |Vup. (3.6) 

uGHl(n) Jb^{xo) Jo. 


We observe that if r ^ 0 then obviously e(r,a:o,r) ^ 0. Thus we use the 
following notation. 

Definition 3.6. We denote 


r(xo, r) := sup {r G R : e(T, xo,r) ^0} . 

Now we discuss some basic properties of the quantities that we have just 
defined. 

Lemma 3.7. The quantity introduced in Deftnition AS. dl is finite, namely 

z[xe),r) G [0,+oo). 

Proof. Let cf G C^{Br) with ||(()||l 2 (b^) = 1, and let u{x) := (f){x — Iq). 
Then |l'«||L 2 (n) = ||m||l2(b^(3,j,)) = = 1, and 

e(T, xo,r)i^T — / |Vup > 0 

J Br{xo) Jq J Br 

provided that r > /g |V(()p. □ 

Lemma 3.8. For any ti ^ T 2 we have that 


e{T2,xo,r) - e(Ti,xo,r) G [0,r2 - n]. 

Proof. Fix e > 0. For any i G {1,2}, there exists M(i,e) G i7{(0), with 
l|w(i,e)||L 2 (Q) = 1 such that 


e{n,xo,r) < £ + Ti 


I Brixo) 


U 


2 

(i,e) 
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Therefore 


e{T2,XQ,r) - e{Ti,XQ,r) ^ T2 j 

JE 


Br{xo) 




^ n / W(l,s) - 

JbAxo) 

^ -e, 


/ |VM(i,e)P - e(ri,a;o,r) 

Jn 

/ |V-«(i,e)|^ - e(ri,a;o,r) 
Jn 


and 


e(Ti,xo,r) - e(T2,xo,r) ^ n / “(2 e) “ / |V'«(2,£)P - e(r2, xq, r) 

> (n - T 2 ) [ 

Jb, 


“(2, e ) “ ^ 


I Br{xo) 

> -{t2-Ti) / W(2,e)-e 

= -{T2-Ti)-e. 

The desired result now follows by taking e as small as we wish. 
Corollary 3.9. If t I T{xo,r), then e{T,xo,r) —>■ 0. 

Proof. Suppose not, i.e. there exists a sequence 

Tj > T(xo,r) 

with Tj —>• t{xo, r) as j —>■ + 00 , such that 


|e(T,-,xo,r)| > a, 


□ 

(3.7) 

(3.8) 


for some a > 0. We claim that 


e{Tj,xo,r) a. (3.9) 

We prove it by contradiction: if not, by (13.811 . we would have that e(Tj,xo, r) ^ 
—a. Thus, we set 

Ta ■■=T[xo,r) + -. 

We notice that > xfxo^r), therefore, by Definition 13.61 we have that 

e(Ta,xo,r) > 0. 

In addition, we have that > r, if j is large enough, thus we make use of 
Lemma l3.81 and we obtain that, for large j, 

0 + a < e{Ta,xo,r) - e(T,,xo,r) ^ Ta - Tj. 

Taking the limit in j, we conclude that 

/ X ® 

a^Ta- T{xo,r) = -. 
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This is a contradiction and (EH) is established. 

Also, by Definition 13.61 we know that there exists a sequence fj ^ lixoTr) 
with Tj —>• r(a;o,r), such that e{fj,xo,r) < 0. Accordingly, by (j3.9|) . 


e{Tj,xo,r) - e{fj,xo,r) ^ a. (3.10) 

Notice that tj > t{xo , r) ^ fj and 

lim Tj — fj = T{xo,r) — t{xo, r) = 0. 

Thus, by Lemma 13.81 

lim e(Tj,xo,r) — e(fj,xo,r) ^ lim Tj — fj = 0. 

t —>■ + 00 t —> + oo 

This is in contradiction with (|3.10l) and so the desired result is proved. □ 

Before stating and proving the main theorem of this subsection, we investi¬ 
gate the behavior of T{xo,r) under scaling. 

Proposition 3.10. Fix s' € (0,1). There exists a constant r* := T*(s',n) such 
that 

z{xo,r) ^ T»(s', n) 

for every xq € Q and r > 0 such that Br{xo) C 12. 

Proof. We claim that 


/ +^c(sM2)r2*'||Vu||i.(^), (3.11) 

JBrixo) 

for some constant c(s', 12) > 0. Once (13.111) is proved, one can finish the proof of 
the desired result by arguing as follows. One sets t*(s', 12) := l/c(s', 12). Then, 
for every r < t-*(s', 12) (i.e. for every r ^ l/(c(s', 12) r^'' )), one has that 


/|Vu|"-r/ ^ [ \Vu\^ -[ 

Jn JBAxa) Jn c(s',12)?-^® 


> 0, 


Br{xo) 


where the latter inequality is a consequence of the claim (I3.11|) . This gives that 
e{T,xo,r) ^ 0 for any r < r*(s', 12), and so, by Definition 13.61 we have 

that T{xo,r) ^ r*(s',12), thus proving the desired result. 

Due to these observations, it only remains to prove (13.111) . To this scope, we 
observe that, given p > 2, by the Holder inequality with exponents | and 
we have 

/ + < (w„r")V ||u||ip(n). 

JBrixo) 

Therefore, the claim in (13.111) is established if we show that there exists p > 2 
such that ^ 

ll«llL(0) ^ C(s', 12,p) +' II Vu||i.(o), (3.12) 
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for some C{s', V,,p)> 0. So, now it only remains to prove (13.121) . To this goal, 
we deal separateljo with the cases n = 2 and n ^ 3. 

We start with n ^ 3. In this case, we denote by p := > 2 the Sobolev 

conjugate exponent of 2. Notice that = 2 and the Sobolev inequality 

(see e.g. formula (7.26) in [53]) bounds ||u||^p(q) with (7(11) ||Vu||^ 2 (q), for 
some (7(n) > 0. Hence, if we denote by Dq > 0 the diameter of H, we have that 

ll^llLp(a) = ^^ll^llip(n) ^ Dq IIVu|||2(q), 

and estimate (j3.12l) follows in this case. 

For the case n = 2, we observe that 

lim -—- = 1 > s', 

p->+oo p 

so we can choose an even integer p = p(s') G (2, +oo) large enough such that 

P-2 


> s . 


P 


Also, the critical Sobolev embedding (see e.g. formula (7.38) in 


exp 


i{x)\_y 


,ci||Vu||i2(Q)y 

for suitable ci, C 2 > 0. Then, since 


dx < C 2 |H|, 


(3.13) 
) yields that 

(3.14) 


+°° +k 






we deduce from (|3.14p that 

for some C{ft,p) > 0. Therefore 

ll^llij>(n) ^ C"(f2,p) ||Vu|||2(n), 

for some C'{il,p) >0. As a consequence, if Dq > 0 is the diameter of il, 


(p-2)r. 


r " l|M||LP(n) 

^C'{n,p)Do 
This completes the proof of (13.121) when n = 2. 


2(^-0 2.' 


llVull 


L2(n) ■ 


□ 


®The case n 3 is simpler because the Sobolev conjugated exponent 2* = 2n/(n — 2) is 
not critical. Indeed, in this case the parameter s' does not play much role. 


25 





















Theorem 3.11. Let r, t > 0. Consider the family of distributions Ct = 
TXBr(xo) O'lT'd a corresponding family of stationary solutions Ur € Hq{LI), that is 

—Aut = (ctt- — Ur)uT ■ 

U T i L{xo,r), then Ur 0 uniformly. 

Proof. First of all, we notice that 

Ur ^ T, (3.15) 

thanks to Lemma [2.131 Now we fix £ € (0,1) and we claim that 

l|Wr||L3(0) ^ £, (3.16) 

provided that r is close enough to T(xo,r). To establish this, we test the equation 
against Ur itself, and we obtain that 

f |VMt-P = [ {(Tr - Ur)ul = T ( ul- f ul, 

JQ JQ J Brixo) Jq 

which in turn gives 

l|wr||i3(o) = ul = T ul- < e(T,xo,r), 

Jci j Br{xQ) Jq, 

thanks to (13.61) . This and Corollary 13.91 imply (I3.16p . 

Now we set g{x) := {dr — Ut)ut. Notice that —Aur = g in 12 and, by (I3.15|) 
and Lemma 13.71 

\g\ < (o-r + Ur) Ur ^ 2 t Ur ^ 2(r(xo, r) + 1) Ur Cq Ur, 

for some Co > 0 independent of r, as long as r is sufficiently close to T(a;o,r). 
In particular, by (13.151) and (13.161) . 

ll5llL"+3(n) ^ Co < Cl < Ci£^, (3.17) 

for some Ci > 0. Moreover, using the Holder inequality with exponents 3 
and 3/2, 

2 

ll“^llL3(n) = ^ ^ (/j I|wr|li3(n), 

therefore, recalling (13.161) and (I3.17L 

l|Wr||L2(0) + ||g||Ln+3(n) ^ se + Ci£"+3 ^ C2£"+3 , 

for some C 2 > 0. We combine this information with Theorem 8.15 of [23] (used 
here with / := 0 and q := 2(n + 3) > n), thus we obtain that 

||'ST||L°°(n) ^ C (||MT||L2(n) + ||5lU3(n)) ^ CC 2 , 

for some C > 0, as long as r is sufficiently close to T(a:o, r), which is the desired 
claim. □ 
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Corollary 3.12. Fix s' G (s, 1). Let r,T>0. Assume that 


r < 


^ Ct{s,Bi)T4s'M )\^^ 


(3.18) 


where C'u(s, i?i) is the Poincare-Sobolev constant in (12.41) and T*(s',n) is given 
by Provosition \3.10[ 

Consider the family of distributions = TXBrCo)- Then there exists r > 
T(xo,r) such that both the reverse Poincare-Sobolev condition in (1^ and the 
mismatch condition in en) are satisfied. 

Proof. By taking r large enough, one can easily fulfill (12.81) . This and Lenima [275] 
guarantee the reverse Poincare-Sobolev condition in (12.61) . 

In particular, by Theoreni l2.10l we can consider the solution Ur correspond¬ 
ing to the resource CTt-. 

Now we fix 

£e('o 


r. O! , ^ • (3-19) 

2j.2(s'-s) I ^ ^ 

Thanks to Theorem 13.Ill we can choose r sufficiently close to :r(a;o, r) such that 
||■u^-||L~(n) ^ r~'^‘^e. Therefore, for every x G Br(xo), we have that 

ar(x) — Ur{x) ^ (Jr{x) — r~'^‘'E > T{xo,r) — r“^®e. 

From this and Proposition 13.101 we have that, for every x G Br{xo), 
arix) — Ur(x) ^ ^(s', LI) — r~^‘'E . 

So, recalling (13.191) . 


inf (ar(x) — Ur(x)) > 
xGx^Brixo) 


r r*(s', LI) 


Thus, from Remark 12.31 and (13.181) . we obtain 

1 1 j,-2s' ^2(s'-s) 


C^{s,Brixo)) r2®C'ii(s, Bi) Cii(s,Bi) 


r-^^' C^{s,Bi)T,(s',Lt) 


< 


2Q(s,Bi 


r-2«'r,(s',F!) 


< inf {aT.{x) — Ur{x)) 

x€.Br{XQ) 


This establishes the mismatch condition in (13.11) . 


□ 


From Proposition l3.2l and Corollarv l3.121 it follows that we have constructed 
another example for which the equilibrium (mt-jO) is linearly unstable, confirm¬ 
ing again Theorem 11.11 Once again, this example corresponds to a resource 
that is unevenly spread in the environment, and the nonlocal diffusion may 
compensate such unbalanced distribution of resource. 

As a final observation, we would like to stress that most of the techniques 
discussed in this paper are of quite general nature and can be efficiently exploited 
in similar problems with different species and different dispersive properties. 
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4 A purely nonlocal phenomenon 


Goal of this section is to rely the stability of stationary points of type (0, v) with 
Theorem ll.2l and show how, with our arguments, there is no hope to prove an 
analogue of Theorem 11.11 for {0,v). We include the proof of Theorem 11.21 and 
clarify that it is a purely nonlocal feature. 

The linearization of the system (USD at (0,{i) gives 


L(o,i){u,v) 


A + (cr — {;) 0 

—V —(—A)® + (cr — 2v) 



Thus, any instability result would be a consequence of an inequality of type 


15(0,U) (’^*7 0) 



v)u\ > 0 . 


This means that, if we want to run the same argument that we did in Section 
[3] for (m, 0), then we have to hnd an analogue for the mismatch condition (13.11) . 
Roughly speaking, we need to know that, at least in certain circumstances, 
the amount of leftovers of the dominant population v exceeds a given constant, 
depending on the size of the domain. But, around a stationary point of type 
(0, {;), the nonlocal population v tends to exhaust all the available resource a in 
the domain O. This claim is motivated by Theorem 11.21 because formula (11.71) 
states that the population Ug locally fits with any given resource, up to an 
arbitrarily small error estimated by dn}. Of course we are neglecting the 
Dirichlet boundary condition on 

We are now left with the proof of Theorem 11.21 

Proof of Theorem \l.‘A By Theorem 1.1 in [20) . we know that we can approx¬ 
imate cr by a s-harmonic function in Bi\ namely, we have that there ex¬ 
ist > 1 and Ue £ C^{Bi) (~l C®(]R"') satisfying p.8l) . (II.9p and 

lk-Me||c'=(Si) ^ (4-1) 

Now we define 

(Te := Ue- (4.2) 

In this framework, formula (11.61) follows from (14.11) and (14.21) . Moreover, by (II.8p 
and (14.21) . 

{(Je{x) - Ue{x)) Ue(x) = 0 = {-AYu^ix), 

for any x G Bi, which proves (11.101) . □ 

We stress that Theorem ll.2l is only due to the nonlocal feature of the equation 
and it does not have any local counterpart, as pointed out by the next result. 

Proposition 4.1. Let M > 0. Let a € C^(Bi) with 

a(x) ^ M for any x G Bi/ie 
and a(x) ^ 1 for any x G Bi \ i?i/io ■ 
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Then, there exists Mq > 0 and e > 0 such that, for any M ^ Mq, ifue S C^(Bi) 
satisfies 

Ik- o-£||c 2 (Bi) ^ ^ k-3) 

and Ue G C^{Bi) satisfies 

- Aue(x) = (fTek) - Ue{x)) Ue(x) for any X G Bi, (4.4) 


then 


\\ue - CTeWc^B^) > £ (4-5) 

In particular, the local counterpart of Theorem \1.2\ is false. 

Proof. Suppose by contradiction that for every e > 0 there exist and Ug 
satisfying not only (14.31) and (14.41) . but also 

Ike - cre||c2(Bi) ^ £. 

From (14.31) and (14.51) . we know that 

Ike - 0 -|koo(Bj) ^ ||ug - cr||c 2 (Si) < 2£. (4.6) 


As a consequence, 


lkelk“(Si) ^ 2 + |kllc 2 (Bi) < e'er, 

for some Ca > 0, possibly depending on the fixed resource cr. This, (H3]) 
and (14.41) give that, in B\, 

I Azigi ^ |(Tg zig I keI ^ £. 

Thus, the weak Harnack inequality (see e.g. Theorem 8.18 in [53]) gives that 

lk£lki(Bi/ 4 ) ^ ei ( inf Ue + Cae], (4.7) 

for some constant Ci > 0. Now, by (jUj) and dSI), we see that zzg(x) ^ M — 2e 
in BijiQ and therefore 


lkelUi(Si/ 4 ) ^ / Ue{x)dx^ C 2 {M-2e), (4.8) 

for some constant C 2 > 0. Similarly, from (jUj) and (ir31) . we have that ztg ^ 
1 + 26 : in Bi \ i?i/io and therefore 

inf Ug ^ 1 + 2s < 2. (4.9) 

-Bi/s 

By inserting (14.8|) and (14.91) into (14.71) we obtain that 

M-2e^ Cs (2 + e). 


for some C 3 > 0. Thus, we take M ^ Mq := SC's. This fixes a and gives that 

C's ^ dVi — 2 C 3 ^ 2e + Ca (2 + e) — 2 C^ = (2 + C^Ca) s. 

By taking e small, we obtain a contradiction and we complete the proof of 
Proposition 14.11 □ 
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5 Further comments on stability and nontrivial 
solutions 


Of course, the results presented in this paper do not aim to exhaust the variety 
of scenarios offered by the analysis of local and nonlocal competing species. 
In particular, further investigations about existence and local/global stability 
of equilibrium solution are desirable, also with the aim of establishing under 
which conditions local and nonlocal strategies are convenient for the evolution. 

In particular, while we focused here on the local stability (i.e. whether or not 
a mutation of strategy turns out to be persistent for small times), the strategic 
question for biological population in competition for large times is mostly related 
to global stability. 

The question of attractors for the global dynamics is related to the regular¬ 
ity properties of the semiflow and to the associated maximum and comparison 
principles. For this reason, though not directly used in this paper, we present 
here in detail a general comparison principle for a single fractional equation 
(the general case of systems deserves a separate analysis, also due to the lack of 
cooperativeness between biological species, see e.g. formula (7) in [^, and we 
plan future further investigation along the lines of [25l [28] ): 

Lemma 5.1. Let T > 0 and consider a locally Lipschitz function f. Let v and 
w be bounded and continuous solutions of 

dtv + {—/S.yv + f{v) > dtw + {—AYw + f{w) (5.1) 


on M" X (0, T], with v{x, f) ^ w{x, t) for any x £ M" \ and any t £ [0, T\, and 
v{x, 0) ^ w{x, 0) for any x £ fl. 

Then v{x, t) ^ w{x^ t) for any for any x £ K" and any t £ [0, T], 

Proof. The proof is of classical flavor, see e.g. Proposition A.5 in [Hj. By 
possibly iterating the argument, it is enough to prove the result up to a small 
time, hence, without loss of generality we may assume that 




1 

4(M-f 1)’ 


(5.2) 


where M ^ 0 is the (local) Lipschitz constant of / ~ more precisely, we take M 
such that 


f{w{x,t) -ri)- f{w{x,t)) ^ M| 77 | for any rj £ [- 1 , 1 ]. 


(5.3) 


Now we suppose, by contradiction, that the claim were false. Then, it would 
exist (x,f) £ n X [0,r] such that v{x,t) < w{x,t). We define 

w{x, t) — v{x, t) T 
4 ’ I 


e := mm 


(5.4) 


and W := V — w + et + e^. Notice that 

IF(x, 0) = v{x, 0) — w{x, 0) -b > 0 


and 


W{x,f) = v{x, t) — w{x, t) + et + e'^ ^ v{x, f) — w{x, t) + 2e < 0, 
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thanks to (lOl) . Hence, there exists 2* := G H x [0, t] C H x [0,T] such 

that W{x,t) > 0 for any x G H and any t G with = 0. 

In particular, IH(a;,t*) ^ 0 = and so the integrodifferential defi¬ 

nition of the fractional Laplacian gives that (—A)®lH(a;*, t*) ^ 0. 

Also, lH(a;*,t) ^ 0 = lH(a;*,t*), and thus i9tlH(a;,,t*) ^ 0. In addition, we 
have that 

v{z^) = lH(z,) -I- u>(z,) — et* — = w{z^) — ry*. 


where 77, ;= et* -|- G 


0 , 


, thanks to (15.21) and (15.41) . 


2(Af-|-l)_ 

As a consequence of these observations, and recalling (15.11) and (15.3L we find 
that 


0 ^ dtW{z,) + {-Ayw{z,) 

= {ptv + {-AYv) {zY - {ptw - {-AYw) {zY} + e 

> f{w[zY) - f{v[zY) + £ 

= f{w{zY)-f{w{zY-r]Y+£ 

^ -M'q^: + £ 

> 0 , 

which is a contradicition. □ 

In addition, we remark that the theory developed in the previous pages also 
allows us to investigate the stability of nonlocal species. For instance, one sees 
that large resources allow both local and nonlocal populations to stem from 
the pure equilibria. More precisely, if a is larger than the first classical and 
fractional eigenvalue, then the reverse Poincare-Sobolev condition in (12.61) (and 
its fractional counterpart (I2.7|) l are satisfied. 

So, in case system (HH) possesses two unstable pure equilibria (fi, 0) and (0, w), 
a positive mixed equilibrium (tt*(x), ?;*(x)) may arise. 

Of course, a detailed analysis of all these circumstances in a general setting 
and a careful check of the nontrivial details involved by the dynamics associated 
to the flow go beyond the scope of this paper, but we refer also to [7] for a series 
of examples which carefully compare local and nonlocal behaviors of biological 
populations in terms of the size of the domain and of the sparseness of the 
resources. 
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